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We review the mechanism of spin-lattice coupling in relieving the geometrical 
frustration of pyrochlorc antiferromagnets, in particular spinel oxides. The 
tetrahedral unit, which is the building block of the pyrochlore lattice, under- 
goes a spin-driven Jahn- Teller instability when lattice degrees of freedom are 
coupled to the antiferromagnetism. By restricting our considerations to dis- 
tortions which preserve the translational symmetries of the lattice, we present 
a general theory of the collective spin-Jahn- Teller effect in the pyrochlore lat- 
tice. One of the predicted lattice distortions breaks the inversion symmetry 
and gives rise to a chiral pyrochlore lattice, in which frustrated bonds form 
helices with a definite handedness. The chirality is transferred to the spin 
system through spin-orbit coupling, resulting in a long-period spiral state, as 
observed in spinel CdC^O^. We discuss explicit models of spin-lattice cou- 
pling using local phonon modes, and their applications in other frustrated 
magnets. 



1.1 Introduction 

As explained in the introductory chapter by Chalker, sufficiently strong frus- 
tration in a magnet results in a large degeneracy of its ground-state man- 
ifold. Prime examples of this behavior are Heisenberg antiferromagnets on 
the kagome [1,2] and pyrochlore [3,4] lattices with interactions restricted to 
nearest-neighbor sites. In the classical limit of a large spin S, the ground 
states of these magnets exhibit very high, continuous degeneracies and pos- 
sess numerous zero modes, which correspond to moving the system through 
its manifold of ground states [5] . The pyrochlore antiferromagnet represents 
a particularly striking example of high ground-state degeneracy: at least half 
of its spin- wave modes have zero frequencies in any collinear ground state [4] . 

A large degeneracy means enhanced sensitivity to perturbations, even 
when these are nominally weak. In this chapter we will consider a coupling 



2 



O. Tchornyshyov, G.-W. Chern 



between spins and the underlying lattice, which has its origin in the depen- 
dence of the exchange integrals on the atomic positions, J(r!,r 2 ) Si • S2, and 
is known as magnetoelastic exchange [6]. In the pyrochlore antiferromagnet, 
this coupling lifts the degeneracy of the classical ground states and induces a 
symmetry-lowering distortion of the lattice, in analogy with the spin-Peierls 
effect in antiferromagnetic spin chains [7] . A spin-Peierls-like phase transition 
has been observed in several antiferromagnetic spinels where the magnetic 
ions form the pyrochlore lattice [8-10]. 

The problem of coupled spins and lattice degrees of freedom in a py- 
rochlore antiferromagnet is reminiscent of the collective Jahn- Teller effect [11] 
in crystalline solids. We therefore begin the discussion by studying the Jahn- 
Teller distortion in a tctrahcdral "molecule" with four spins, which is the 
structural unit of the pyrochlore lattice. A symmetry-based analysis will be 
supplemented by models with specific spin-phonon coupling mechanisms. We 
will then extend the analysis to an infinite lattice, to examine some of the 
possible ground states of the classical spin system. In concluding this chapter, 
we will test the theory of spin-phonon coupling on the example of CdCr2C>4, 
a frustrated Heisenberg antiferromagnet with S = 3/2 spins residing on the 
pyrochlore lattice. 

1.2 Spin-driven Jahn- Teller effect in a tetrahedron 



Fig. 1.1. Representation of four spins in the corners of a regular tetrahedron. 

Consider the basic structural unit of the pyrochlore lattice, a cluster in 
the shape of a regular tetrahedron with four spins of length S at the corners 
(Fig. 1.1). The Heisenberg exchange energy depends on the total spin of the 
cluster, S to t = Sj=i S», according to 
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For antiferromagnetic exchange (J > 0), this energy is minimized when the 
total spin of the cluster is 0. The degeneracy of the ground state is thus equal 
to the number of distinct spin-singlet states in a system of four spins. 

For spins of length 5, the number of linearly independent singlet ground 
states is 25 + 1. Indeed, the pair of spins 1 and 2 can have a total combined 
spin S12 ranging from to 25. The same is true of the spin pair 3 and 4. A 
state with a total spin 5tot = can be obtained by combining states with 
5i2 = 5,34 = 0, 5i2 = 534 = 1, and so on, up to 5i2 = 534 = 25. Thus 
one observes that, for any spin length 5, the ground state of the four spins is 
degenerate. The high symmetry of the cluster means that, in accordance with 
the Jahn- Teller theorem, 3 the ground-state energy can be lowered through a 
distortion. 

Assume for simplicity (a more general case will be considered below) that 
the exchange coupling J between spins i and j has a dependence on their 
separation, = |r, — r ^ | , given by 



where is a reference distance. At zcroth order in the displacements Sri, we 
recover the unperturbed Heisenberg Hamiltonian (1.1) with 25+1 degenerate 
singlet ground states. 
The first-order term, 



lifts the degeneracy of the ground-state manifold. As long as the displacements 
involved remain small enough to satisfy J 'Sr^ <C J, excited states with 5 to t > 
at energies J and higher may be neglected. Thus it is necessary to determine 
the energy levels of H\ in the Hilbert space of the singlet ground states. 

1.2.1 Generalized coordinates and forces 

The perturbation Hamiltonian (1.3) depends on the atomic displacements 5vi, 
which comprise 4x3=12 degrees of freedom. However, not all of these will 
influence the exchange energy of the spins. As one example, this does not 
change under rigid translations of the tetrahedron (3 modes) or under global 
rotations (a further 3 modes). The remaining 6 modes represent various defor- 
mations of the four-site cluster. To facilitate further analysis, we classify these 
modes in terms of the irreducible representations (irreps) of the tetrahedral 
point group T d [13]. 

The 6 modes belong to three irreps of T d . The breathing mode (irrep A) 
leaves the symmetry of the tetrahedron fully intact. A doublet of tetragonal 



J(r y ) = Jin,) + J'if^Snj + J" (f + ... 



(1.2) 




(1.3) 



i<j 



3 It is important to note that the degeneracy is not caused by the symmetry of 
time reversal, so the conditions of the theorem are fulfilled [12]. 
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and orthorhombic distortions, Q E — (Qf,Qf), transforms under irrep E. 
Finally, a triplet Q T2 = (Q^ 2 , Q^ 2 , Q J 2 ) , transforming as irrep T 2 , elongates 
and compresses opposing bonds; equal-amplitude superpositions of the triplet 
components yield trigonal distortions. The coordinates of these modes can be 
expressed in terms of Cartesian displacements of the spins with the coefficients 
listed in Table 1.1. 



Table 1.1. Coefficients relating the 6 distortions of a tetrahedron to the displace- 
ments <5r; of its vertices. The reference frame is shown in Fig. 1.1. 
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Table 1.2. Coefficients relating the bond elongations Srij to the distortion coordi- 
nates Q A , Q E = (Qf,Qf), and Q T2 = (Q? , , Q^). 
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By expressing the changes in bond lengths, Srij, in terms of the generalized 
coordinates, we reduce the perturbation Hamiltonian (1.3) to the form 

a n 

where the index a enumerates the irreps and n its components. The variable 
f" is the generalized force which is conjugate to the coordinate and has 
the same symmetry properties. J a/ is a coupling constant appropriate for the 
irrep a. 

For illustration, the breathing mode Qa couples to the spin operator 
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which is invariant under all symmetry operations of T^. Further, up to a trivial 
factor, this is the unperturbed Hamiltonian (1.1) and so has the same value 
in any of the degenerate ground states. Consequently, the term —J A Q A f A 
in the perturbation Hamiltonian (1.3) produces a trivial energy shift of the 
degenerate ground states, but does not split them. 

The triplet mode and the associated triplet force also do not induce a 
splitting. To demonstrate this, we note that Qi 2 couples to the operator 

fl 2 = (S 2 ■ S 3 - Si ■ S 4 )/A (1.6) 

which vanishes in any singlet state. 4 In the presence of an applied magnetic 
field, the triplet forces cannot be neglected because of the nonzero total spin 
iStot- The triplet forces and the corresponding trigonal distortions play an 
important role in the stabilization of the half-magnetization plateaus observed 
in some spinel chromitcs (Sec. 1.3.1). 

The only two modes involved in the splitting of the ground state are the 
components of the doublet (Qf , Qf") of tetragonal and orthorhombic distor- 
tions. These couple, respectively, to the spin operators 



ft 



E Si • S4 + S2 • S3 + S2 • S4 + Si • S3 — 2Si • S2 — 2S3 • S4 



12 

jts 1 ~+ "i S2 "j (17) 



E S2 • S4 + Si • S3 — Si • S4 — S2 • S3 



In what follows, we omit the irrep superscript because only the doublet E 
participates in lifting the degeneracy of the ground-state manifold. 

In addition to the magnetoelastic exchange (1.3), which is linear in the 
distortion amplitude, it is necessary to consider also the elastic energy of the 
deformation. We thus obtain the spin-distortion Hamiltonian 

H = J'(Q • f) + fc|Q| 2 /2 ee J'(Qi/i + Q2/2) + k{Q\ + Q 2 2 )/2, (1.8) 

where k is the elastic constant of the doublet distortion. Having established 
this, the next task is to minimize the energy (1.8) with respect to both the 
coordinates and the spins. 



1.2.2 Four S — 1/2 spins on a tetrahedron 

The problem of four S = 1/2 spins on a deformable tetrahedron was first 
analyzed by Yamashita and Ueda [14]. The ground state of the unperturbed 

4 Indeed, from Si + S 2 + S 3 + S 4 = it follows that (Si + S4) 2 = (S 2 + S3) 2 and 
thus Si ■ S 4 = S 2 • S 3 . 
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exchange Hamiltonian (1.1) is two-fold degenerate. As a basis in this Hilbcrt 
space one may use singlet states with a well-defined total spin on bonds 12 
and 34, S12 — S34 = a, where a — or 1 [15]. In this basis, the force operators 
/1 and /2 are proportional to the Pauli matrices o\ and 03, respectively, so 
that the Hamiltonian (1.8) reduces to 

H = J'(QiCT! + Q 2( 7 3 )V3/2 + k(Ql + Ql)/2. (1.9) 

For a given distortion Q, the ground-state manifold is split into the two 
energy levels 

£1,2 = ±\J'\QV3/2 + kQ 2 /2, (1.10) 

and the energy of the system is minimized when Q = \J'\y/3/(2k). Note that 
this quantity depends on the magnitude of the distortion, Q = \/Q\ + Qf , 
but not on its "direction:" it can be tetragonal, purely orthorhombic, or any 
combination of these. This degeneracy is associated with a continuous sym- 
metry of the Hamiltonian (1.9) that involves simultaneous "rotations" of both 
the distortion coordinates and the Pauli matrices, 

Qi\ ( cos 9 sin#\ f Qi\ ( a\ \ ( cos 9 sin#\ / a\ 
Q 2 )^ \-sin9 cos6» ) \Q 2 ) ' \a 3 ) ^ ^-sin^ cos<9 ) \a z 

(1.11) 

The invariance of the energy under this transformation does not reflect an 
underlying symmetry, and applies at the level of approximation used here. 
Terms of higher order in Q break this symmetry to leave only a three-fold 
degeneracy, as may be expected on symmetry grounds [14]. The lowest-order 
anharmonic term allowed by the symmetry is proportional to 

Q x Q v Qz = (~\Qi + ^Q^j (-\Qi - ^-Q*J Qi = jQ 3 cos3a, (1.12) 

where Q Xl Q y , and Q z measure tetragonal distortions along the respective 
axes and a is the polar angle in the (Qi,Q2) plane. Depending on the sign 
of the cubic term, it favors three distinct ground states at a = 0,±27r/3 
or at a — 7r,±7r/3. The former "vacua" have spin singlets on two opposing 
bonds (e.g. S12 — S34 = for a — 0) while the latter have spin triplets 
on two opposing bonds (S12 — S34 = 1 for a = w). These ground states 
exhibit valence-bond order, which violates the point-group symmetry of the 
cluster but not the SO (3) symmetry of the exchange interaction. The two- 
component order parameter (1.7), introduced by Harris, Berlinsky and Bruder 
[16], measures the differences in spin correlations on the different bonds. 

The ground states of the cluster exhibit a tetragonal lattice distortion 
along one of the three major axes with Q = — J'(f)\/3/(2fc). If J' < 0, as 
would be expected for direct antiferromagnetic exchange, the tetrahedron is 
flattened (elongated) in a ground state with triplets (singlets) on opposing 
bonds. 
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For spins of length S > 1/2, the analysis proceeds by a similar route [15]. 
The lowest-order perturbation (1.3) yields three degenerate singlet ground 
states with the largest possible spins on two opposing bonds, such as 5 i2 = 
£34 = 25, and a tetragonal distortion (a flattening of the tetrahedron for 
J' < 0). This result is most easily understood in the classical limit S — > 00, 
to which we turn next. 

1.2.3 Four classical spins on a tetrahedron 

For classical spins, the total energy E(f , Q) (1.8) can be minimized in two 
steps. We minimize it first with respect to the distortion Q at a fixed f. 5 A 
minimum is achieved when Q = — J'f/fc, yielding the energy 

25(f) = - J' 2 / 2 /(2fc), (1.13) 

whence the total energy is minimized by states with the largest magnitude of 
the force doublet f . Thus it is necessary to quantify the magnetoelastic forces 
in those states of the ground-state manifold with S^ot = 0. 

The domain of attainable f values forms a regular triangle in the (/1, / 2 ) 
plane (Fig. 1.2). Its three corners correspond to the three distinct collincar 
states with four satisfied bonds (Si • Sj = —S 2 ) and two frustrated ones 
(S,-Sj = +S 2 ). States elsewhere on the perimeter of the triangle have coplanar 
spins. 

Not unexpectedly, the doublet force is maximized (and the total energy 
minimized) in the collinear states. Indeed, in such states antiparallel spins 
attract and parallel spins repel each other with forces F = — J' Si - Sj of the 
largest possible magnitude, | J' IS* 2 . The large forces result in large distortions, 
yielding a large decrease in the total energy. Thus in the classical limit one 
expects collinear ground states in which the tetrahedron is flattened along one 
of the principal axes for J' < 0. 6 

1.2.4 Color notation and other useful analogies 

We find it convenient here to introduce an analogy with the color triangle, 
where the vertices correspond to the primary colors red, green, and blue, the 
midpoints of the edges to the secondary colors cyan, magenta, and yellow, 
and the center to the absence of color. If we color bonds perpendicular to the 
major axes a, 6, and c respectively red, green, and blue, then the color of the 
state in Fig. 1.2 reflects the color of the frustrated bonds. 

The collinear ground states obtained for classical spins provide a simple 
rationalization for the valence-bond states found for quantum spins of length 

5 This method cannot be applied to quantum spins because the operators /1 and 
fi do not commute [15], and so their values cannot be measured simultaneously. 

6 If J' > 0, the spins are still collinear but the tetrahedron is elongated along the 
same axis. 
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Fig. 1.2. The domain of attainable values of the force doublet f = (/i, $2) (1-7) for 
classical spins. Dashed lines indicate frustrated bonds (Si • Sj > 0). Reprinted with 
permission from Ref. [15]. 

S > 1/2 [15]. Quantum states maximizing the spins of opposing pairs (for ex- 
ample S12 = S34 = 2S) are the analogs of the collinear classical configurations 
(Si = S2 = —S3 = —S4). 

Lastly, we recall that the spin-distortion Hamiltonian (1.8) was derived for 
a simplified model of exchange in which the Heisenberg interaction strength J 
depends only on the separation of the two spins. This approximation is good 
for direct spin exchange [17], which is the dominant exchange interaction in 
the chromium spinels ZnCr204 [18] and CdCr204 [19], as well as in some 
other chromium antiferromagnets [20]. However, in other situations J may 
exhibit a more complex dependence on atomic displacements, such as the 
very sensitive bond-angle-dependence of superexchange. Fortunately, the form 
of the spin-distortion coupling derived above (1.8) is robust, as can be seen 
by symmetry considerations: group theory guarantees that there are no other 
invariant terms of the same order in f and Q. In the general case, J' represents 
a linear combination of exchange derivatives. 
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1.2.5 Spin-Jahn- Teller effect on a triangle 

Another well-known lattice producing strong frustration in an antiferromagnet 
is the kagome geometry [1,2], a network of corner-sharing triangles in a plane, 
and its three-dimensional variant, the hyper kagome lattice [21]. It is natural to 
ask whether spin-lattice coupling is also an effective mechanism for relieving 
frustration in this case. We answer this question by considering the building 
block of the kagome lattice, a triangle with three spins. 

Classical Heisenberg spins interacting through antiferromagnetic Heisen- 
berg exchange minimize their energy by making angles of 120° with one an- 
other. An analysis along the lines of Sec. 1.2 shows that the correction to the 
exchange energy from the magnetoelastic term is quadratic in the spin dis- 
placements: the linear term cancels. Without this linear term, a spontaneous 
distortion does not occur. The absence of the linear term can be understood 
simply from the magnetoelastic forces between the spins: the three forces, be- 
ing proportional to the scalar products of the spins, are equal in a state where 
the spins make equal angles with each other. These symmetrical forces only 
shrink the triangle without distorting it. 

The argument against the Jahn- Teller distortion fails if the quadratic term 
in the magnetoelastic correction is negative and large enough to overcome the 
purely elastic cost of the distortion, a scenario proposed recently by Wang 
and Vishwanath [22] . In our view, however, empirical evidence indicates that 
a Jahn- Teller instability of this sort would be a rare exception. The strength 
of the quadratic magnetoelastic term can be estimated from the splitting 
of degenerate phonons in antiferromagnets with a spin-induced Jahn- Teller 
distortion [18]. Such splittings, observed in ZnCr204 [18], CdCr204 [19], and 
MnO [23], do not exceed 12% of the phonon frequencies, which indicates the 
dominance of the purely elastic contribution. 



Fig. 1.3. Valence-bond ground states of three spins S = 1/2 with antiferromagnetic 
Heisenberg exchange interactions. The ellipses designate singlet bonds. 

While the spin- Jahn- Teller distortion of a triangle appears unlikely in the 
classical limit, the opposite extreme — quantum spins of length S = 1/2 — 
are a completely different situation. The ground state of three such spins in- 
teracting through antiferromagnetic Heisenberg exchange is any state with a 
total spin Sa = 1/2. Such a state is fourfold degenerate: part of this degener- 
acy is of the Kramers type, as the projection of the total spin on an arbitrary 
axis can be = ±1/2, and there is an additional twofold degeneracy related 
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to the symmetry of the triangle. This degeneracy can be understood in terms 
of valence-bond states, in which two of the spins on the triangle form a singlet 
bond (a quantum dimer) while the third remains free. Figure 1.3 shows three 
such states, although they are not mutually orthogonal; in fact only two of 
these states are linearly independent. 

The presence of a non-Kramers degeneracy leads to the spin- Jahn- Teller 
effect on a triangle with spins S = 1/2. The analysis is similar to that for 
four S = 1/2 spins on a tetrahedron (Sec. 1.2.2), with three distinct ground 
states. Depending on the sign of the cubic term, the distorted triangle will 
have either two long bonds and one short bond with a singlet on it, or two 
short bonds and one long bond with a triplet. 

1.3 Models with local phonon modes 

The symmetry-based analysis of the previous section is basically exact regard- 
less of the underlying microscopic model for the phonons. In this section we 
review some specific models of spin-lattice coupling based on local phonon 
modes and also discuss their applications. 

Probably the simplest situation is the "bond-phonon model," in which the 
exchange integral (1.3) and also the elastic energy depend only on the bond 
length Tij . The elastic energy is approximated by the sum of individual terms 
k5rfj/2. After integrating out the bond variables Srij, the model generates a 
biquadratic spin Hamiltonian 



which clearly favors collinear spin configurations. Because of its simplicity, 
this model has been applied in numerous studies of frustrated magnets. As 
one example, Bccca et al. found using the bond-phonon model that magne- 
toelastic coupling leads to a spin-Pcierls transition in the frustrated J\-Ji 
antiferromagnet on the square lattice [24,25]. This study may be relevant 
to the nature of the transition to a phase of collinear order observed in the 
quasi-2D antiferromagnet I^VOSiC^. 

On the pyrochlore lattice, the bond-phonon model gives a ground state 
with 3 w -fold degeneracy, where N is the total number of tetrahedra: each 
tetrahedron can be flattened along one of the 3 major axes, independently 
of the other tetrahedra. A more realistic phonon model can be formulated in 
terms of the independent displacements of each atom, with the bond lengths 
determined from these atomic displacements by Sr^ = (iij — Uj) • fy. This is 
known as the site-phonon (or Einstein phonon) model of spin-lattice coupling 
[26,27], in the simplest version of which the elastic energy is approximated by 
a sum of individual terms fc|ui| 2 /2, an assumption which leads to a constant 
dispersion similar to the long- wavelength limit of optical phonons. In addition 




(1.14) 



(ij) 
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to Eq. (1.14), after integrating out the displacements the model introduces 
couplings between bond variables, 

-(J ,2 /2k) ]T (f« • f ifc ) (Si ■ Sj) (Si ■ S fc ). (1-15) 

Because of these coupling terms, coherent long-range distortions are possible 
in this model. Using the site-phonon model, Wang and Vishwanath [22] showed 
that a zigzag collinear order of the triangular antifcrromagnet CuFeC>2 is 
a ground state of the spin-lattice Hamiltonian. However, because the Fe 3+ 
ion in this compound has spin S = 5/2 and L = 0, resulting in a rather 
small magnetic anisotropy, a relatively large spin-lattice coupling is required 
to induce the zigzag collinear order from the non-collinear 120° ground state 
of Heisenberg spins on this lattice. 

In the context of the pyrochlore lattice, the interaction term of Eq. (1.15) 
corresponds to an antiferromagnetic coupling between the force doublets of 
nearest-neighbor tetrahedra, KJ2{ a p)^a • fp, where the coupling constant 
K > 0. As a result, neighboring tetrahedra tend to be flattened along dif- 
ferent directions (e.g. (100) and (010)). While this reduces the total number 
of ground states, it still leaves a very high accidental degeneracy which grows 
exponentially with N. 

1.3.1 Half-magnetization plateau in ACr 2 4 spinels 

Local phonon models provide an explicit description with which one may 
study the spin-lattice stabilization of the half-magnetization plateau observed 
in some spinel chromites. The low-temperature magnetization curves of the 
spinels CdCr 2 4 and HgCr 2 4 exhibit a sharp transition into a wide plateau 
region where the magnetization is equal to half its saturation value [28,29]. 
Each tetrahedron in the half-plateau state is in one of the four 3-up-l-down 
collinear spin configurations. While thermal and quantum fluctuations act 
in general to favor collinear spins, and indeed in some cases help to sta- 
bilize the magnetization plateaus in frustrated magnets, the observed half- 
magnetization plateau in spinels arises most likely due to a coupling between 
the spins and the lattice, a scenario also supported by the observation of a 
spontaneous structural distortion accompanying the half-plateau transition. 

Because the total spin S t ot is nonzero in the presence of a magnetic field, 
coupling of the moments to the singlet (A) and triplet (T 2 ) phonon modes 
can no longer be neglected. Still, when the applied magnetic field is weak, 
the distortion is such that the crystal retains a tetragonal symmetry. The 
spins develop a canted antiferromagnetic order with two frustrated bonds 
and four antiferromagnetic bonds. As the field strength increases, the doublet 
force If^l = (4/V3) (S 2 — S 2 ot /16) decreases as a result of the increasing 
total spin. At a critical field, the trigonal distortion, accompanying a 3-up-l- 
down collinear spin configuration maximizing the triplet forces, fj 2 = \/2S 2 , 
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becomes more favorable energetically. The tetragonal distortion thus gives 
way to the trigonal one through a discontinuous transition. 

Using the bond-phonon model, Penc and co-workers [30] obtained a classi- 
cal phase diagram for the spin- Jahn- Teller effect in a single tetrahedron. The 
transition to the collinear 3:1 states takes place at a field strength H w 3 J and 
J' 2 /k > 0.05 J. A general symmetry analysis also showed that the collinear 
3:1 states are always stabilized over a finite range of magnetic field, provided 
that J E ' < 2J T ? [30]. 

Similar to the ground-state manifold at zero field, the model with indepen- 
dent bonds retains an extensively degenerate manifold of half-magnetization 
states, because the trigonal distortion of individual tetrahedra can be along 
any of the four (111) axes. This accidental degeneracy is lifted by the ad- 
ditional coupling term (1.15) introduced by the site-phonon model, which 
favors an antidistortive coupling between neighboring tetrahedra. Using this 
rule, Bergman et al. [27] showed that a unique spin configuration (up to dis- 
crete lattice symmetries) with a quadrupled, cubic unit cell is the ground 
state of the system in the half-magnetization plateau. The resulting space 
group, P4 3 32, is consistent with the X-ray diffraction pattern of the spinel 
HgCr 2 4 [31]. 

1.4 Collective Spin- Jahn- Teller effect on the pyrochlore 
lattice 

Attempts to extend the analysis of the spin- Jahn- Teller effect on a single 
tetrahedron to an infinite pyrochlore lattice encounter a conceptual problem: 
there are infinitely many phonon modes coupled to the spins (one may ex- 
pect two for every tetrahedron). There are also technical difficulties: detailed 
knowledge of the crystal's elastic properties is required. As a result of these 
difficulties, the problem lacks a general solution. 

Some progress may, however, be made through the use of local phonon 
models, as described in the previous section. Still, a massive accidental de- 
generacy remains in the ground states of these simplified models. Further, the 
lattice modes in real crystals are plane waves, and thus a lattice distortion 
involves only a small number of phonons with specific lattice momenta. For 
example, the distortion in ZnCr204 shows superlattice peaks in a diffraction 
experiment with wavenumbers [\\\] [10]. This does make it possible to take 
an alternative, phenomcnological approach in which only a small number of 
lattice modes is considered. Such an approach was taken by Tchernyshyov 
et al. [15], who considered the simplest case where spin displacements pre- 
serve the translational symmetries of the lattice and break only point-group 
symmetries. 

The pyrochlore lattice is made up of tetrahedra of two different orienta- 
tions. Because all tetrahedra of the same orientation are related by lattice 
translation (which is assumed to remain a good symmetry), it is necessary 
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to consider only two tetrahedra of opposite orientations A and B (Fig. 1.5). 
The symmetry group must be expanded from by including the inversion 
I exchanging the two sublattices of tetrahedra, T d ® I — Oh [13]. The irreps 
remain largely unaltered, with the exception of a newly added parity index, 
which enters because these are either even (g) or odd (u) under the inversion. 

At linear order in the displacements, the only modes which couple to the 
spins are the doublets E g and E u . The former represents an overall tetragonal 
or orthorhombic distortion of the lattice, while the latter is an optical phonon 
with wavenumber q = that distorts tetrahedra of opposite orientations in 
exactly opposite ways {e.g. by flattening tetrahedra A and elongating tetra- 
hedra B by the same amount and in the same direction). These modes can 
be expressed in terms of linear combinations of distortions on tetrahedra of 
types A and B, 

Q A + Q B Q A ~Q B 

Q 3 = - , Q" = - J 4 (1-16) 

The spin-lattice energy (1.8) generalizes to 

E(i A , t B ,Q A , Q B ) = ,r(Q A -{ A + Q B -t B ) + k g \Q3\ 2 /2 + k u \Q u \ 2 /2, (1.17) 

where k g and k u are the elastic constants of the even and odd E doublets. 
Minimization with respect to the lattice modes Q 9 and Q u yields the energy 
as a function of the spin variables in the form 



E(i A ,i B ) 



K a \f A + i B \ 2 K u \i A -f B \l 

(1.18) 



4 4 
(K g + K u )(\i A \ 2 + \i B \ 2 ) 



4 

(K g -K u )({ A -f B ) 
2 

where we have introduced the effective magnetoelastic exchange couplings 
K — J' 2 Ik 

Jv g.u u I ™g,u • 

The second line in Eq. (1.18) is the result familiar from Eq. (1.13): the 
magnitude of the doublet force f is maximized on tetrahedra of both sublat- 
tices. Thus one expects to find states with collinear spins and (for J' < 0) 
tetrahedra flattened along one of the three (100) directions. The third line in 
Eq. (1.18) represents a coupling between the f variables of the two sublattices, 
whose consequences depend on which of the two lattice modes is softer. 

If K g > K u (Eg mode softer), the energy is minimized when f A and f B 
are in the same corner of the force triangle (Fig. 1.2). Tetrahedra of both 
sublattices are flattened along the same (100) direction, so that only an E g 
distortion is present. The spin configuration is shown in Fig. 1.4. The magnetic 
unit cell coincides with the structural one. Because we are considering an 
0(3)-symmetric Heisenberg model, the global orientation of the spins can be 
arbitrary 
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Fig. 1.4. Magnetic order in a state with a softer E g mode. Frustrated bonds are 
shown as dashed lines. The lattice is flattened uniformly (for J' < 0). Reprinted 
with permission from Ref. [15]. 

If, on the other hand, K g < K u (softer E u mode), the two f vectors are 
located in different corners of the triangle, so that tetrahedra of types A and 
B are flattened along two different (100) directions, giving six possible ground 
states. The distortion is a superposition of the E u and E g modes. The presence 
of the even mode is understood in a straightforward manner: if tetrahedra of 
type A are flattened along (100) and tetrahedra of type B along (010), the 
lattice is on average elongated along (001). The presence of the E u component 
of the distortion means that the inversion symmetry is broken spontaneously. 
This has important consequences for the magnetic order, shown in Fig. 1.5, 
as we discuss in detail below. Here we note only that frustrated bonds form 
left-handed spirals in this particular state. 

1.5 Collective Jahn- Teller effect in CdCr 2 04 

The normal spinels ACr2 04 (where A is a nonmagnetic Zn, Mg, Cd, or 
Hg ion) are strongly frustrated antiferromagnets exhibiting the spin-induced 
Jahn- Teller distortion. The magnetic Cr 3+ ions forming the pyrochlore lat- 
tice have electronic configuration 3d 3 . The oxygen octahedron surrounding 
the chromium ion splits the 3d levels into a high-energy e g doublet and a 
low-energy ti g triplet; the former are unoccupied while the latter are singly 
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occupied, with the 3 electrons having parallel spins. Thus the orbital degree of 
freedom is quenched and the spins form a nearly isotropic magnetic moment 
with total spin 5 = 3/2 and a gyromagnetic ratio of g = 2 [8]. Interactions 
between the spins are mediated mostly by direct antiferromagnetic exchange 
between neighboring Cr sites [18] (with the exception of HgC^CU, where this 
contribution is comparable to the ferromagnetic superexchange term mediated 
by the oxygen ions). 

All of these compounds have a strongly correlated paramagnetic state well 
below their Curie- Weiss temperatures, 0, and order magnetically at a tem- 
perature Tjy -C 0. The phase transition is discontinuous and is accompanied 
by a lattice distortion [8-10]. The lack of orbital degrees of freedom points to 
a magnetoelastic origin for the lattice distortion. 

CdCr 2 4 is the best-understood system to date. In the distorted state 
below TV it shows no superlattice peaks [9] , indicating that the translational 
symmetry of the high-temperature cubic phase (space group Fd3m) remains 
unchanged. The point-group symmetry is lowered: the lattice exhibits a tetrag- 
onal distortion with lattice constants a = b < c (an overall elongation). The 
low-temperature structure was identified by Chung et al. as the pyrochlore 
lattice with a uniform tetragonal elongation (space group IAi/amd) [9], i.e. a 
pure E g distortion. However, as listed next, there are good reasons to believe 
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that this distortion also involves the staggered component E u , which breaks 
the inversion symmetry and lowers the space group down to IAi 22. 

1. The dominance of direct antifcrromagnetic exchange between adjacent 
chromium spins means that the exchange constant decreases with increasing 
ionic separation, i.e. J' < 0. If the distortion were of pure E g type, the 
crystal would flatten along one axis, yielding a = b > c in contradiction to 
the experimental data [9]. As explained in the previous section, an _E u -driven 
distortion would lead to an overall elongation of the lattice, a = b < c, in 
agreement with the experiment. 

2. An E u distortion breaks the inversion symmetry of the crystal, mak- 
ing the crystalline lattice chiral. Indeed, the elongated (frustrated) bonds in 
Fig. 1.5 form spirals of one helicity. Spin-orbit coupling in the form of the 
Dzyaloshinskii-Moriya interaction would then spread the chirality from the 
lattice to the spins, generating a spiral magnetic order. The observed mag- 
netic order in CdCr 2 04 is in fact a weakly incommensurate spiral [9]. 

1.5.1 Spiral magnetic order in CdCr 2 4 

Chung et al. reported an incommensurate magnetic order with magnetic Bragg 
peaks at q = 2-7r(0, S, 1) in a crystal with an elongated c axis (a = b = 
0.995c) and 6 = 0.09. The magnetization lies in the ac-plane. Because the 
incommensurability 6 is small, the magnetic order can be understood as a 
commensurate state with q = 27r(0, 0, 1) twisting slowly along the ^-direction. 

The same authors found two such structures which would be consistent 
with the data they obtained from elastic neutron scattering. One of the pro- 
posed ordering patterns is derived from the commensurate state shown in 
Fig. 1.5, and is precisely what one expects when the magnetoelastic effect is 
driven by the E u phonon. The other candidate is derived from an "orthogo- 
nal" state where the spins on every tetrahedron are oriented, for example, in 
directions +x, — x, +y, and — y. Such a state is very hard to obtain through 
the spin-driven Jahn- Teller effect [15], and no other justification for this state 
is currently known. 

The small value of 5 makes it possible to treat this incommensurability 
as the result of a weak perturbation on top of the Heisenberg exchange and 
magnetoelastic coupling. Chern et al. [32] derived possible magnetic spiral 
states that arise when the Dzyaloshinskii-Moriya interaction is added to these 
two energy terms. These authors found two candidate solutions, one of which 
is entirely consistent with the data of Chung et al. This analysis and its results 
are described in the next section. 

1.5.2 Theory of spiral magnetic order 

The Dzyaloshinskii-Moriya (DM) interaction gives a contribution 

m 
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to the Hamiltonian, where the sum extends over pairs of nearest neighbors. 
This term is allowed on the ideal pyrochlorc lattice, where the bonds are not 
centrosymmetric, the precondition for a non- vanishing coupling constant Dy. 
At the same time, the high symmetry of the lattice constrains the direction 
of this vector [33]: for a bond oriented along the [110] lattice direction, the 
vector must point along [110] (Dy = (±D,^fD, 0)/\/2)- The value of Dy on 
any other bond is then found through the symmetry transformations of the 
system. In a collinear magnetic state, the expectation value of -Hdm is zero, 
but its contribution can be lowered by twisting the spins into a spiral. 

The pitch of the spiral is determined by the competition between the DM 
coupling strength D and a spin stiffness. In most antiferromagnets, the spin 
stiffness is set by the exchange interaction, but the pyrochlore antiferromagnet 
with nearest-neighbor exchange is an exception: the large degeneracy of its 
ground state leads to a vanishing stiffness. As a result, the ground state in 
the presence of the DM interaction is not a long-range spiral but is rather a 
commensurate state with noncollinear spins [33]. 

The presence of magnetoelastic interactions changes the situation. Recall 
that the spin-induced Jahn- Teller effect selects a state with collinear spins. 
This tendency towards collinear states results in a finite spin stiffness. As a 
result, the pitch of the spiral is a quantity of order D/ K, where K = J' 2 /k is 
the effective magnetoelastic interaction. 

A problem where the nearest-neighbor exchange J, the magnetoelastic 
energy scale K, and the DM coupling D are each arbitrary is hard to solve 
analytically. However, it simplifies if these energy scales are well separated, 
the conventional hierarchy being 

JS 2 > KS i > DS 2 . (1.20) 

A quantitative analysis indicates that in CdCr 2 C>4 these scales are of similar 
magnitude (as discussed at the end of this section), but still have the expected 
order, JS 2 > KS 4 > DS 2 . Thus while a theory based on the assumption 
(1.20) may not provide a quantitative account of magnetic order in CdCr 2 04, 
it presents at minimum a good point of departure for understanding it. 

Minimization of the dominant term, the exchange energy (1.1), yields a 
constraint that the total spin be zero on every tetrahedron. The remaining 
degrees of freedom of a single tetrahedron can be parametrized using the 
staggered magnetizations 

_ Si — S 2 — S3 + S4 _ —Si + S2 — S3 + S4 _ —Si — S2 + S3 + S4 
Ll 4S ' L2 4S ' Ls " 4S 

Because each spin belongs to two tetrahedra, the staggered magnetizations on 
one sublattice of tetrahedra determine completely those of the other sublat- 
tice. We use the staggered magnetizations of sublattice A, {L^ 4 }, to express 
the magnetizations of sublattice B, and, except in cases of possible confusion, 
suppress the sublattice index to simplify the notation. 
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Even for a single tetrahedron, the three staggered magnetizations are not 
completely independent. Vanishing of the total spin in the ground state makes 
them mutually orthogonal and imposes on their lengths the constraint [34] 

3 

U=Li%, 1,1, 5 {j , ^i? = l. (1.21) 

i=l 

The lengths Li parametrize the angles between the spins and are related to 
the bond doublet f by 

.h = 2S 2 (Ll + Li-2L 2 3 )/V3, h = 2S 2 {Li — L\). (1.22) 

Thus five parameters arc required to describe the magnetic state of a tetra- 
hedron in the ground state of Eq. (1.1): three Euler angles for the triad {lj} 
and two further parameters for the bond doublet, e.g. f = (/i, / 2 ). 

The magnetic order of CdCr 2 04 in the commensurate limit S — > 
(Fig. 1.5) has staggered site magnetizations: L 2 = L 3 = and Li = e* q r ni, 
where q = 27r(0, 0, 1) and ni is an arbitrary unit vector. In terms of the three 
staggered magnetizations, the DM term for a single tetrahedron is 

Ebm = -DS 2 (a • L 2 x L 3 + b • L 3 x Li + c • Li x L 2 ). (1.23) 

It is easy to see that the DM energy is exactly zero for the collinear state. For 
cither sign of D, this term can be made negative by allowing a small compo- 
nent of L 2 or L 3 , which describes a twisting spin configuration. The lowering 
of the DM energy (1.23) is accompanied by an increase of the magnetoelastic 
energy. However, further analysis shows that the former is linear in L 2 and L3 
while the latter is quadratic, so that such a twisting distortion always occurs. 

To pass to a continuum description, we express the rapidly oscillating unit 
vectors li = e* q r in terms of a slowly varying orthonormal triad {n^}, and 
use the length constraint (1.21) to eliminate L\ by 

Li = - h±h? j ftl e iq-r ) L2 = L2 A2 e iq-p ) Ls = ^ n 3 e iq r . (1.24) 

At this point, the magnetic structure is described in terms of five slowly 
varying fields, L 2 , L 3 , and the triad {ni}. 

The number of independent fields can be further reduced by examining 
tetrahedra of sublattice B: the vector of total magnetization on each B tetra- 
hedron must vanish, giving three more constraints. The total spin of a B 
tetrahedron centered at r B = r A + (1/4, 1/4, 1/4) is given by 

M B (r B ) = Sx(r A + &1 ) + S 2 (r^ + a 2 ) + S 3 (r A + a 3 ) + S 4 (r A ), 

where a x = (0,1/2,1/2), a 2 = (1/2,0,1/2), and a 3 = (1/2,1/2,0) are primi- 
tive lattice vectors (the centers of the tetrahedra form a diamond lattice; the 
Bravais lattice is fee). Expressing the spins {S^} in terms of the staggered 
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magnetizations {L,} and using the gradient expansion, one obtains, to the 
lowest orders in L 2 , L 3 , and the gradients, the constraint 

M* = L 3 n 3 -i|± = 0. (1.25) 

Thus it is clear that L3 and A3 are determined by the gradient Vni. 

In a similar way, expressing the staggered magnetizations on sublattice B 
to lowest order in Vni leads to 

L 1 -L 2 n 2 - i — , L 2 -n 1; L 3 - — — . (1.26) 

Substituting these expressions into Eq. (1.23) and adding contributions from 
tetrahedra of both types yields 

£ DM = - ^n, • (a x — + b x — - c x — ). (1.27) 

The DM energy contains terms linear in gradients of hi, indicating that the 
spins are unstable against the formation of spiral configurations. As an ex- 
ample, the first term ni • a x dhi/dx favors a magnetic order with the unit 
vector ni perpendicular to the a-axis and spiraling about it. 

As discussed previously, the spin stiffness has its origin in the magnetoelas- 
tic energy. For simplicity, here we consider only distortions due to E u phonons 
and neglect the effect of E g distortions (a procedure definitely appropriate in 
the limit K g <C K u ). The linear decrease of the DM energy due to gradients 
of fii must be balanced by the increase in magnetoelastic energy, which on 
symmetry grounds must be quadratic in Vni. From Eq. (1.18), the increase 
of magnetoelastic energy is E mc = —K u u ■ Su/2, where K u = J' 2 /k u and 
u = f A — i B . The unperturbed odd doublet is Uo = 45* 2 (0,1). Because we are 
describing the spiral magnetic order in terms of the staggered magnetizations, 
it is convenient to use Eq. (1.22) for the calculation of 5f A and 5f B . Retaining 
terms to second order in Vni leads to 



^ 1 \ 2 + f^ 1 \ 2 +2(^ 1 V L " + 4L 2 

dx ) \ dy J \ dz ) 2 2 dz 2 



.(1.28) 



Because the DM energy (1.27) does not depend on L 2 , minimization of the 
total energy with respect to L 2 affects only the magnetoelastic term (1.28) 
and yields 

thus L 2 is also eliminated. The minimized magnetoelastic energy is 
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The total energy of the spiral state, now expressed as a functional of the 
vector fields rii(r) and n 2 (r), is the sum of Eqs. (1.27) and (1.30). Its min- 
imization yields a second-order partial differential equation. While we are 
unable to find the most general solution to this equation, we can find three 
highly symmetrical spiral solutions in which the spins remain perpendicular 
to, while twisting about, one of the (100) axes. As one example, a spiral state 
along the b axis is described by ni = (cos 9(y), 0, sin 9{y)) and has total energy 

E = -(DS 2 6' +K u S 4 e' 2 )/4, (1.31) 
where 0' = d6/dy. Minimization of this quantity gives the pitch of the spiral, 

e '= 2 « 5 =2£jp- (L32) 

Equation (1.29) implies that L^ 1 = 0; the A tetrahedra therefore have coplanar 
spins spanned by two orthogonal Neel vectors, and L^. The angles between 
spin pairs are given by 614 = 623 = 2L 3 . From Eq. (1.25), this angle is related 
to the pitch by 2L 3 = irS. The spiral magnetic state has the structure 



Lf = cos(7r5/2) (cos(27r<5y), 0, sin(27rfo/)) e 

■ A 
J 2 



\>2 = 0, (1.33) 



2niz 



L3 4 = sin(7T(5/2) (-sin (271%), 0, cos {2n8y)) e 

producing Bragg scattering at wavevector q = 27r(0, ±5, 1), while the ordered 
magnetic moments lie in the ac-plane [Fig. 1.6(b)]. All of this is consistent 
with the experimental data of Ref. [9] . 

It is worth noting that the distorted crystal structure preserves certain 
symmetries on interchanging the A and B sublattices, such as inversion at 
a Cr site followed by a 7r/2 rotation in the a6-plane. However, the magnetic 
order described here breaks these symmetries. From Eq. (1.26) one has Lf = 
Lf = 0, which means that every tetrahedron on sublattice B has collinear 
spins, whereas the spins of the A tetrahedra are twisted into a (weakly) non- 
collinear state. This disparity between the two sublattices should result in 
different distortions of the two types of tetrahedra, thus further lowering the 
symmetry of the crystal. However, the magnitude of the additional distortion 
is expected to be small because the degree of non-collinearity is small, 5 <C 1. 

Spiral states in which the spins rotate about the a-axis [Fig. 1.6(a)] can 
be obtained similarly, by using the ansatz fix = (0, cosd(x), sm9(x)). The 
resulting solution can also be obtained from Eq. (1.33) through symmetry 
operations which exchange the two sublattices of tetrahedra, such as inversion 
on a Cr site. This spiral produces a magnetic Bragg peak at q — (=p<J, 0, 1). 

Finally, there is a third spiral solution, shown in Fig. 1.6(c), where the 
spins twist about the c-axis. In this state, which is not related to the other 
two solutions by any symmetry, the magnetic Bragg peak occurs at q = 
27r(0, 0, 1 + S) and both sublattices have tetrahedra with coplanar, rather than 
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Fig. 1.6. Three symmetrical spiral magnetic structures minimizing the energy. The 
spins are perpendicular to and twist about the a-axis (a), 6-axis (b), or c-axis (c). 
Dashed lines indicate frustrated bonds. The crystal is viewed along a (100) direction. 
Reprinted with permission from Ref. [32]. 



collinear, spins. That this spiral state has the same energy as the previous two 
is a coincidence: its total energy is also given by Eq. (1.31) with 0' = d6/dz. 
This degeneracy is lifted when other perturbations, such as further-neighbor 
interactions, are taken into account. CdCr204 has a significant third-neighbor 
antiferromagnetic exchange interaction, which acts to favor strongly the states 
with spirals twisting along the a- or 6-axis [32] . 

In closing this section, we comment on the assumption (1.20) of well- 
separated energy scales in the problem. Ab initio calculations [32] yield a 
nearest-neighbor exchange energy JS 2 = 1.1 meV and a magnetoelastic en- 
ergy K U S 4 — 0.76 meV. The strength of the DM interaction can be estimated 
from the measured pitch of the spiral, S = 0.09, using Eq. (1.32), which gives 
DS 2 = 0.21 meV. While the three energy scales are not vastly different, they 
do appear in the correct order of decreasing magnitude, JS 2 > K U S 4 > DS 2 . 



1.6 Summary and open questions 

The spinel compound CdCr204 provides an opportunity to test our under- 
standing of the ground state of the Heisenberg antiferromagnet on the py- 
rochlore lattice. When the lattice degrees of freedom are included, both the 
selected magnetic order and the lattice distortion are in agreement with a 
theoretical model [15] based on two vibrational doublets of the crystal, the 
q = optical phonon E u and the uniform lattice distortion E g . The model 
ties the incommensurate nature of the spiral magnetic order to a spontaneous 
breaking of the inversion symmetry in the crystal, which has not yet been 
observed directly in CdCr 2 04. 

The magnetoelastic phase transition between the high-T correlated para- 
magnet [35] and the low-T ordered phase remains poorly understood. It is 
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strongly discontinuous in both ZnCr 2 04 [8] and CdCr 2 04 [9], with both the 
lattice distortion and the ordered moment reaching their T = values im- 
mediately below the ordering temperature. A Landau free-energy approach 
based on the spin-Peierls order parameter (1.7) appears to be the only can- 
didate approach available for modeling the underlying physics [15]. However, 
it is not evident that this phenomenology can provide a full description. In 
CdC^Cu, the order parameter has the E u symmetry and does not allow for a 
cubic invariant in the free energy, thus excluding the most obvious cause for a 
discontinuous phase transition. The phase transition may yet turn out to be 
first-order if the even (i? 9 -symmetric) order parameter is nearly soft [15], but 
this question has not yet been clarified. In fact it remains unclear whether the 
valence-bond variables (1.7) represent a good choice of the order parameter 
for these systems: the low-T phase in both ZnCr 2 04 and CdCr 2 04 is mag- 
netically ordered, suggesting that a type of staggered magnetization may be 
a more appropriate choice. 

Finally, a realistic model of this phase transition must take into account 
the entropy of the correlated paramagnetic state [4] . The very high entropy of 
the disordered phase may be responsible for the discontinuous nature of the 
phase transition, as has been demonstrated for the case of lattice models with 
large numbers of flavors [36]. 
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